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Abstract. We analyse pre-Main Sequence evolutionary tracks for low mass stars with masses m < 1.4 Mq based 
on the Baraffe et al. (1998) input physics. We also extend the recent Chabrier et al. (2000) evolutionary models 
based on dusty atmosphere to young brown dwarfs down to one mass of Jupiter. We analyse current theoretical 
uncertainties due to molecular line lists, convection and initial conditions. Simple tests on initial conditions show 
the high uncertainties of models at ages <J 1 Myr. We find a significant sensitivity of atmosphere profiles to the 
treatment of convection at low gravity and T e g < 4000 K, whereas it vanishes as gravity increases. This effect adds 
another source of uncertainty on evolutionary tracks at very early phases. We show that at low surface gravity 
(log g <C 3.5,) the common picture of vertical Hayashi lines with constant T e s is oversimplified. The effect of a 
variation of initial deuterium abundance is studied. We compare our models with evolutionary tracks available in 
the literature and discuss the main differences. We finally analyse to which extent current observations of young 
systems provide a good test for pre-Main Sequence tracks. 
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1. Introduction 

The development of a new generation of stellar evolu- 
tion models based on the accurate coupling between in- 
terior and atmosphere models yielded a major advance 
in the description of very-low-mass stars (VLMS; Baraffe 
et al. 1995, 1997, 1998, BCAH98) and substellar objects 
(brown dwarfs BD and giant planets GP; Burrows et al. 
1997 for objects with T cff < 2000K; Chabrier et al. 2000a, 
CBAH00). One of the main advantages of such models 
over previous generation models is the direct compari- 
son between theory and observation in colour-colour and 
colour-magnitude diagrams (CMD). Several observational 
tests now assess the validity of the theory devoted to 
the description of low-mass (<1M Q ) astrophysical ob- 
jects (see Allard et al. 1997; Chabrier and Baraffe 2000, 
for recent reviews). Although some discrepancies between 
models and observations still remain, uncertainties due to 
the input physics are now significantly reduced. 

In previous papers, we have developed evolutionary 
models of VLMS and BDs (BCAH98; CBAHOO), taking 
into account the most recent improvements of the physics 
describing the interior (equation of state for dense plas- 
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mas, screening factors) and the atmosphere (molecular 
opacities, formation of dust) devoted to the analysis of 
objects with an age t> 100 Myr. Comparison between 
observations and models for very young objects (t < 100 
Myr) are more uncertain, in particular from the observa- 
tional viewpoint, mainly for two reasons : (1) extinction 
due to the surrounding dust modifies both the intrinsic 
magnitude and the colours of the objects and (2) spectra 
of very young objects (t < 1 Myr) may still be affected by 
the presence of an accretion disk or circumstellar mate- 
rial residual from the protostellar stage. On the theoret- 
ical side, an important source of uncertainty at the very 
early stages of evolution is the choice of the initial condi- 
tions, resulting from prior protostellar collapse and accre- 
tion phases. The simple picture of non-accreting objects 
contracting from large initial radii, as used in our models, 
is clearly an idealized description of reality (Stahler, Shu 
& Taam 1980; Stahler 1983, 1988; Hartmann, Cassen & 
Kenyon 1997 and references therein). 

In spite of these uncertainties, numerous surveys are 
devoted to very young clusters because of the potential 
detection of substellar objects down to planetary masses 
(Zapatero et al. 2000; Lucas et al. 2001), and the deter- 
mination of the very (sub)stellar initial mass functions 
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(IMF), free of dynamical evolution effects that affect older 
clusters. Consequently, a wealth of data for low mass ob- 
jects with ages spanning ~ 1-10 Myr is available and pro- 
vides the basis for a better understanding of the early 
phases of stellar evolution. Given the reliability of the 
present theory for VLMS and BD, a confrontation of the 
models with the observational data of such very young 
objects provides important insight into the early epoch of 
protostellar collapse and accretion phases. For this pur- 
pose, however, a correct appreciation of the model uncer- 
tainties and limitations is required. 

For this reason, this paper is specifically devoted to the 
analysis of evolutionary models by BCAH98 and CBAH00 
at early phases of evolution (t < 100 Myr). The input 
physics is briefly recalled in section 2, followed by a dis- 
cussion of uncertainties specific to evolutionary tracks at 
early ages (§3). Our models are compared to other evo- 
lutionary tracks available in the literature in §4 and the 
confrontation to observations is presented in §5. 

2. Model description 

2.1. Evolutionary tracks 

The models analysed in the present paper are based on the 
input physics already described in BCAH98 and CBAH00. 
Both sets of models use the same ingredients describ- 
ing the stellar interior but use different sets of atmo- 
sphere models, which provide the outer boundary condi- 
tions and the synthetic spectra. The BCAH98 evolution- 
ary tracks are based on the non-grey atmosphere mod- 
els by Hauschildt, Allard and Baron (1999a). These mod- 
els are dust-free and are appropriate to the description 
of objects with effective temperatures T e g > 2300 K. The 
CBAH00 models are based on atmospheres including the 
formation and opacity of dust (Allard et al. 2001, here- 
after DUSTY models ). As illustrated in CBAH00, dust 
must be taken into account in order to explain the near- 
IR colors of late M-dwarfs and L-dwarfs. The latter mod- 
els are thus more appropriate to the description of ob- 
jects with T off <; 2300 K. As emphasized in CBAH00, the 
DUSTY models are not appropriate for the description 
of spectral and photometric properties of methane dwarfs 
(T c ff < 1600K), which require a different treatment of dust 
(the so-called COND models in CBAH00 and Allard et al. 
2001). Evolutionary models for such cool objects will be 
presented elsewhere (Chabrier et al. 2001, in preparation). 

The BCAH98 grid covered a mass range from 0.02 Mq 
to 1.4 Mq, for ages > 1 Myr up to the Main Sequence for 
stars. Originally, the CBAH00 grid covered masses from 
0.01 M Q to 0.1 M Q for ages > 10 Myr. In the present 
work, we extend it down to 1 Mj (1O _3 M0) and ages > 
1 Myr. Figure |l| presents the complete grid of models in 
a Hertzsprung-Russell diagram (HRD) from 0.001 Mq to 
1.4 Mq. The small variation of radius with mass and age 
in the substellar regime (see e.g. Chabrier & Baraffe, 2000) 
yields the merging towards very similar tracks below the 
hydrogen-burning limit. Figure || displays the time evo- 



lution of the effective temperature and luminosity for se- 
lected masses. Objects below 2 Mj evolve essentially with 
T c ff < 1600K (see Figs, [j] and 2), even at very early ages. 
Their atmospheric properties are thus better described by 
the COND models. 

The initial conditions of the models are described in 
the next section. The stellar/substellar transition is lo- 
cated at m ~ 0.075 Mq. Below this limit, objects become 
partially degenerate and their nuclear energy produc- 
tion cannot compensate the energy lost by radiation (e.g 
-^nuc < L), which is required to reach the Main Sequence. 
The deuterium burning minimum mass is m ~ 0.012 Mq 
(Saumon et al. 1996; Chabrier et al. 2000b). The initial 
D-burning phase lasts less than 1 Myr for m >0.2M Q , be- 
tween 1 and 5 Myr for 0.05<m<0.2M Q and almost 20 
Myr for a 0.02 M brown dwarf (see Chabrier et al. 2000b, 
for details). 

As shown in CBAH00, grains have little effects on 
the evolution of Lit) and T e s(t), because of the reduced 

~ 1 /3 

dependence of evolution upon opacity, Lit) cx n R , 

T eS (t) cx k~ 1/W (Burrows & Liebert, 1993). We verified 
that the difference between the different TiO and H2O 
molecular line lists used in BCAH98 and CBAH00 mod- 
els (see §3 in CBAH00), respectively, affect essentially the 
outer atmospheric layers, and thus the synthetic spectra 
and colors, but not the deeper atmospheric layers, and 
thus the outer boundary conditions. The effect of these dif- 
ferent molecular linelists on the evolution of the effective 
temperature T c ff (t) and the bolometric luminosity L(t) is 
small, less than 100 K in T e ff and 10% in L at a given age. 
This is illustrated in Figure 2, where tracks based on the 
BCAH98 input physics (solid lines) are compared to the 
CBAH00 tracks (dashed lines) for the same masses (0.075, 
0.02 and 0.012 Mq). As stressed in CBAH00 and Baraffe 
et al. 2001, the computation of more reliable H2O and, to 
a lesser extend, TiO linelists is badly needed to solve this 
shortcoming in the present theory. 

3. Main uncertainties: initial model and 
convection 

3.1. Initial conditions 

As mentioned in the previous section, although shortcom- 
ings still remain in current molecular opacities, the result- 
ing uncertainty on the evolution is small. One of the main 
source of uncertainty for models at early stages of evo- 
lution is the choice of the initial conditions. Most of low 
mass pre-Main Sequence (PMS) models available in the 
literature (D'Antona & Mazzitelli 1994, 1997; Burrows et 
al. 1997; BCAH98; Siess et al. 2000) start from arbitrary 
initial conditions, totally independent on the outcome of 
the prior proto-stellar collapse and accretion phases. The 
initial configuration is that of a fully convective object 
starting its contraction along the Hayashi line from ar- 
bitrary large radii. Evolution starts prior to or at cen- 
tral deuterium ignition, with initial central temperature 
~ 5 10 5 K. According to studies of low-mass protostellar 
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Fig. 1. Evolutionary tracks in the Hcrtzsprung-Russell diagram for masses from 1.4 Mq to 0.001 Mq (dashed lines) 
and ages spanning from 1 Myr to the ZAMS (for stars). Several isochrones for 1, 10 and 100 Myr are indicated by 
solid lines from right to left. The location of the ZAMS for stars down to 0.075 Mq is also indicated (left solid line). 



collapse and accretion phases, such initial conditions are 
oversimplified, and low mass objects should rather form 
with relatively small radii (Hartmann et al. 1997, and 
references therein). Based on spherical accretion proto- 
stellar models, Stahler (1983, 1988) defined a birthline in 
the Hertzsprung-Russell diagram where young objects be- 
come visible. Evolutionary tracks should then start from 
this birthline, which fixes the age t — 0. Ages determined 
from models based on the above-mentioned oversimplified 
conditions should then be corrected accordingly, with sub- 
stantial corrections for systems younger than a few Myr 
(see Palla and Stahler 1999). Furthermore, collapse and 
accretion are unlikely to proceed spherically. Spherical col- 
lapse does not consider angular momentum transport, an 



important issue of the early phases, which affects the sub- 
sequent cooling and formation of the protostar. Hartmann 
et al. (1997) recently illustrate the sensitivity of the birth- 
line locus assuming that accretion proceeds through a disk 
rather than spherically. This work stresses again the high 
uncertainty of assigning ages from HRD positions for the 
youngest objects. 

This analysis demonstrates convincingly that assign- 
ing an age to objects younger than a few Myr is totally 
meaningless when the age is based on models using over- 
simplified initial conditions. As shown in the next section, 
for t < 1 Myr, the evolutionary tracks themselves are sen- 
sitive to the initial conditions, whereas after a few Myr, the 
models converge toward the same track. This is the main 
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Fig. 2. Evolution of luminosity and effective temperature 
as a function of time (in yr) for masses from 1.2 Mq to 
0.002 Mq. Masses (in Mq) are indicated on the right hand 
side of the figure. The solid lines are the BCAH98 models. 
The dashed lines are dusty-models based on the CBAH00 
input physics. 

reason why we provide confidently evolutionary models for 
ages t > 1 Myr, considering that below such ages, models 
are too uncertain. To solve this substantial uncertainty re- 
quires the consistent evolution between the 3D collapse of 
the protostellar phase and the subsequent PMS evolution. 

3.1.1. Effect of the initial radius and the mixing length 
parameter 

In this section we examine the effect of the variation of the 
initial radius (or initial gravity) in the initial conditions. 
For this purpose we have extended our grid of atmosphere 
models to gravities logg < 3.5. We have checked the va- 
lidity of the plane-parallel approximation used in these 
models by comparison with similar atmosphere models 
which take into account the effects of spherical geometry. 
These effects are found to become important only for sur- 
face gravities log g < 2 (Hauschildt et al. 1999b; Allard, 
Hauschildt & Schweitzer 2000), well below the range of 
gravities characteristic of the evolution of young low mass 
objects. 

We first calculate a set of models, labeled (A), with 
initial radii fixed to obtain initial surface gravities log g ~ 



Fig. 3. Effect of the initial radius on the evolution of lu- 
minosity and effective temperature as a function of time 
(in yr) for several masses (indicated near the curves in 
Mq). The solid lines correspond to the first set of mod- 
els (A) with initial gravity log g ~ 3 — 3.5 and the dashed 
lines correspond to models (B) starting with initial gravity 
\ogg ~ 2.5. The mixing length is Z m i X = Hp 

3 — 3.5 and initial thermal time-scales Uh ~ a few Myr. 
For masses above the deuterium burning minimum mass 
(DBMM), to > 0.02 M©, such tracks start at deuterium 
ignition, with central temperature ~ 5 10 5 — 10 6 K, defin- 
ing a Zero-Age-Deuterium Burning Sequence. Such ini- 
tial conditions are similar to that used in BCAH98 and 
CBAH00. A second set of models, labeled (B), starts with 
larger radii such that the initial surface gravity log g ~ 
2.5. These initial models are more luminous than the pre- 
vious ones, with central temperatures below 5 10 5 K and 
initial thermal time-scales t t h ~ 10 5 yr. 

We also analyse the sensitivity of the models to the 
mixing length l m i K , characteristic of the mixing length for- 
malism (MLT) used to describe convection. We used two 
different values of the mixing length parameter amix^] = 1 
and 2 (see §3.2 for justifications). 

The time evolution of L and T c s for models (A) (solid 
lines) and (B) (dashed lines) is displayed in Figure || for 
several masses and for a m ; x = 1. Figure [| shows the im- 
portance of the initial radius on T e ff, and thus on L, during 
the first Myr of evolution. After a few Myr, however, the 

1 a m ix = Imix/Hp, with Hp the pressure scale-height 
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tracks based on different initial radii merge, for a given 
mass. For masses ranging from ~ 0.01 to 0.5 M©, the two 
sets of models start at very different T c g. Models start- 
ing with the lowest gravity are cooler by up to several 
hundreds K compared to initially denser, less luminous 
models with the same mass. Note also that for models 
(B), T c ff increases during the early evolution, under con- 
traction, in contrast to the first set of models with initial 
log g > 3.0. This is the consequence of the different surface 
gravities, which strongly affect the atmosphere profiles for 
T cS = 2200-3500K, as illustrated in Figure |. As shown 
in this figure, for this range of effective temperatures, the 
atmosphere profiles also show an extreme sensitivity to 
the mixing length / m j X at low gravity (log g < 3.5) but 
become rather insensitive to a variation of this parameter 
for log g > 3.5 (see Fig. |b). 

Such a sensitivity of the atmospheric profiles to gravity 
and mixing length in this range of T e g stems from the less 
efficient formation of molecular H2 as gravity decreases. 
Figure || displays the fraction of H2 along the atmosphere 
profiles shown in Fig. [| (below the photosphere). For T e g 
= 2500K and logg = 2.5, the fraction iih 2 of H2 decreases 
rapidly along the inner profile. For a m i X = 1, it becomes 
negligible at optical depth r = 100, where the boundary 
condition for the inner structure is defined (see Chabrier 
& Baraffe 1997, their §2.5). For logg > 3.0 (see Fig. g), 
nn 2 exceeds 20%, even in the deepest atmospheric layers. 

As described in Chabrier & Baraffe (2000, §2.2.2), hy- 
drogen atom recombination yields a rapid increase of the 
H2 Collision Induced Absorption opacity (kcia oc fi| 2 ) 
and a decrease of the adiabatic gradient . Both effects fa- 
vor convection, yielding a flatter atmosphere profile, i.ea 
decrease of T at fixed P. This in turn favors the formation 
of H2 and increases Hh 2 at fixed P, increasing again the 
opacity and the efficiency of convection. This illustrates 
the non-linear response of the atmospheric structure to 
the formation of H2 , as soon as its fraction becomes sig- 
nificant. It explains the huge effect of gravity when this 
latter increases from log g = 2.5 to 3.0 on the atmospheric 
profile, as illustrated in Figure [|a. The decrease of V a d in 
regions of H2 formation favors adiabatic convection, lim- 
iting the extension of super-adiabatic layers and thus the 
sensitivity of the thermal profile to the mixing length pa- 
rameter. This explains the strong dependence of the at- 
mosphere profile on a m i X for log g — 2.5 (see Fig. ||b). At 
higher gravities convection described by the MLT is es- 
sentially adiabatic and almost insensitive to the choice of 
aw (see §3.2). 

A similar behavior is found in the range of effective 
temperatures 2200 K < T cS < 4000K. Above - 4000K, 
the structure is too hot for H2 to form in significant frac- 
tion whereas below T e ff ~ 2200^, the outer layers are 
dense and cool enough for H2 to form efficiently, even at 
log g = 2.5, . 

The drastic modification of the atmospheric profile, for 
<2mix = 1) as gravity increases from log g = 2.5 to log g > 
3.0 explains the different T e fr's for models starting from 
different initial radii (Fig. ||) . For models (B) , starting with 
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Fig. 4. (a) Gravity effect on P - T atmosphere pro- 
files with T e ff = 2500 K and ^ m ; x = Hp. Pressure is in 
dyne.cm~ 2 and temperature in K. The surface gravity 
varies from log g = 2.5 to 4.5, as indicated on the figure. 
The open circles on each curve correspond to an optical 
depth r = 1 . The filled triangles indicate the onset of con- 
vection. The inner atmosphere profiles stop at t — 100, 
where the outer boundary condition for the inner structure 
is defined, (b) Effect of a variation of Z m i x for atmosphere 
profiles with T s = 2500 K and gravity logg — 2.5, 3.5 
and 4.5. Symbols on the dashed curve are the same as in 
(a). 



logg = 2.5, gravity increases as contraction proceeds and 
favors H2 molecular formation, yielding significantly flat- 
ter atmospheric profiles, as mentioned above. This yields 
a significant increase of T c g at fixed mass m as contrac- 
tion proceeds (Chabrier & Baraffe, 2000, §3.2). For a m i X 
= 2, the sensitivity of atmosphere profiles to g is less pro- 
nounced, but still yields up to 200 K differences between 
models (A) and (B) at a given age. 

The effect on evolutionary tracks in a HRD is dis- 
played in Figure [|. Note that evolution along the Hayashi 
line does not necessarily proceed at constant T e ff, be- 
cause of the effects described above. The common picture 
of vertical (constant T e s) Hayashi tracks is therefore an 
oversimplified picture of PMS evolution. As demonstrated 
in the pioneering papers of Hayashi (1961) and Hayashi 
& Nakano (1963), fully convective and adiabatic objects 
which have H2 dissociation zones contract almost verti- 
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3.1.2. Effect of the initial deuterium abundance 
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Fig. 5. (a) Number fraction of molecular H2 as a function 
of P (in dyne.cm~ 2 ) along the profiles displayed in Fig. 
U with T c ff = 2500 K, £ m i x = Hp and different gravities 
log g (indicated near the curves), (b) Same as in (a) for 
'mix = Hp (solid lines) and l m i X — 2Hp (dashed lines). 



cally in the HRD. The very low value of V a d (< 0.1) in 
such zones is responsible for such evolution at constant 
T e g. Taking into account super- adiabaticity in convective 
layers displace the Hayashi line toward lower effective tem- 
perature. This is equivalent to a decrease of a m i x in our 
calculations. Hayashi & Nakano (1963) have also shown 
that omitting the presence of H2 molecules yields an evo- 
lution proceeding at decreasing T c g (from the left to the 
right in a HRD), rather than at constant T c g. Although 
based on a simplified treatment of the atmospheric prop- 
erties, opacities and equation of state, these basic works 
already illustrated the extreme sensitivity of the shape of 
PMS tracks for low mass objects to super-adiabaticity and 
to the degree of H2 formation/dissociation. 

Models (A) and (B) follow the same track for a given 
mass and a m i x , but do not reach the same position at the 
same age. Significant differences appear at ages < 1 Myr 
but vanish after a few Myr. We thus consider 1 Myr as the 
characteristic time required to forget our arbitrary initial 
conditions and below which models are too sensitive to 
input physics and thus too uncertain. 



Deuterium burning plays a key role during the protostellar 
collapse phase (Stahler 1988) and the following ~ 10 Myr 
of evolution. In a previous paper, we focused on the initial 
deuterium burning phase (Chabrier et al. 2000b), adopt- 
ing an initial mass fraction [ 2 Dq] = 2 10~ 5 , characteristic 
of the local interstellar medium (hereafter LISM, Linsky 
1998). In the present section we examine the effect of a 
variation of [ 2 D ] by a factor of 2 on the early stages of 
evolution. Such a variation is motivated by a recent deu- 
terium abundance determination towards quasars, which 
suggests a primordial abundance only slightly larger (by 
less than a factor of 2) than the LISM value (O'Meara et 
al. 2001). 

The main effect of a variation of the initial [ 2 Z?o] is 
the modification of the age at a given L and T c g , with the 
most important effect for masses m < 0.07 Mq (Figure f7j). 
For a standard [ 2 D ] , a 0.07 Mq depletes by a factor of 100 
its initial deuterium content within ~ 3 Myr and a 0.02 
M Q needs ~ 17 Myr (see Chabrier et al. 2000b). If [ 2 D ] is 
increased by a factor of 2, the 0.07 Mq brown dwarf then 
needs ~ 5 Myr and the 0.02 M© requires ~ 26 Myr for the 
same depletion factor. A 0.015 Mq brown dwarf requires 
50 Myr in the standard case and 70 Myr if [ 2 Dq] is twice 
as large to reach the 99% depletion limit. Conversely, if 
[ 2 Dq] is smaller by a factor of 2, the 99% depletion limit 
is reached in 1.7, 12 and 40 Myr for respectively 0.07, 0.02 
and 0.015 Mq. If the initial [ 2 D ] is as small as 2xl0~ 6 , 
e.g a factor of 10 smaller than the standard LISM value, 
all objects with masses m > 0.04 Mq depletes their initial 
deuterium within 2 Myr. 

The effect of a variation of [ 2 £>o] on isochrones in 
CMDs, however, remains small, compared to other sources 
of uncertainties at young ages (extinction, T e g calibration, 
initial models) and can be ignored for this present. If the 
improvement of observable techniques in the future, how- 
ever, allows the detection of deuterium in the atmosphere 
of young objects (Chabrier et al. 2000b), this effect needs 
to be taken into account for a correct age estimate of very 
young clusters. 

3.2. Convection 

One of the major uncertainties affecting the evolution of 
stars with masses to>0.6Mq is due to the treatment of 
convection. These stars show relatively extended super- 
adiabatic outer layers during PMS and MS evolution, of 
which description is extremely sensitive to the adopted 
treatment of convection, i.e., within the MLT formalism, 
to the adopted mixing length (at any gravity). The effect 
of a variation of £ m ; x on evolutionary tracks for solar-type 
stars is well known and is illustrated e.g in Figure 2b of 
Baraffe et al. (2001). Ludwig, Freytag & Steffen (1999) cal- 
ibrated the mixing length parameter a m ; x for these stars 
with 2D hydrodynamical models performed in the param- 
eter space 4300 K < T cS < 7100 K and gravities 2.54 
< log9 ^ 4.74. They found a moderate variation of the 
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Fig. 6. Evolutionary tracks in the Hertzsprung-Russell diagram for masses from 0.2 Mq to 0.01 M , as indicated. 
The solid curves correspond to a m ; x = 1 and the dashed curves to a m [ x = 2. Ages of 0.1 and 1 Myr are indicated by 
respectively crosses and plus for models (A) with initial gravity log g = 3.0-3.5 and by open circles and squares for 
models (B) with initial log g = 2.5. Note that for models (A) the initial position (age=0) is essentially the same as the 
position at 0.1 Myr (cross) and can differ significantly from position of models (B) at t = 0. 



mixing length parameter around typically a =1.5. Note 
that for the present models and input physics ^, l m i x = 
1.9 Hp is the value required to fit the Sun at its present 
age. An increase of a m i x from 1 to 2 yields an increase of 
T e g up to 500K for the highest masses during their PMS 
evolution. 

For masses m <0.6Mq, the extension of the super- 
adiabatic layers retracts appreciably and the transition 
from convective to radiative outer layers is characterized 

2 As stressed in BCAH98, the value of l m i x required to fit the 
Sun depends on the input physics (equation of state, opacities, 
outer boundary condition) 



by an abrupt transition from a fully adiabatic to a radia- 
tive structure with a very small entropy jump. This means 
that during most of the evolution, except at early ages, as 
discussed in §3.1.1 and below, the sensitivity of the evolu- 
tionary models to Z m i X is small. Multi-dimensional hydro- 
dynamical simulations for conditions characteristic of M- 
dwarf atmospheres, T e g < 3000 K, log g — 5, have been re- 
cently conducted by Ludwig, Allard & Hauschildt (2001). 
These simulations confirm the afore-mentioned small en- 
tropy jump found in the 1-D models described by the 
MLT, illustrating the large efficiency of atmospheric con- 
vection for these objects, a direct consequence of the for- 
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Fig. 7. Deuterium depletion curves as a function of age 
in L and T c s for different initial D abundances: [ 2 Dq] — 
10~ 5 (dash-dotted curves), [ 2 D ] = 2xl0~ 5 (solid curves) 
and [ 2 Dq] = 4xl0~ 5 (dashed curve). For each initial value 
of [ 2 -Do]> two curves are shown corresponding to the 50% 
(left curve) and 99% (right curve) depletion limits. The 
symbols on the curves indicate the location of respectively 
0.07 Mq and 0.02 Mq brown dwarfs. 



mation of molecules, as mentioned earlier. Under these cir- 
cumstances, the 3D simulations show that the MLT does 
indeed provide a correct thermal profile, providing a value 
of a m i x > 1, at least for high gravities and relatively old 
objects (t » 10 Myr). 

In §3.1.1, however, we have shown that even below 0.6 
Mq , very young models with gravities log g < 4 can be af- 
fected by a variation of Z m i x . To minimize such uncertain- 
ties, a correct calibration of Z m jx requires the extension 
of the Ludwig et al. (2001) calculations to lower gravi- 
ties. This work is under progress and represents the most 
promising method for an accurate description of convec- 
tion in optically-thin media, through an accurate calibra- 
tion of Z m i x . In contrast, given all the already mentioned 
uncertainties inherent to either observation or models for 
very young objects, a calibration of Z m i X based on a com- 
parison of PMS tracks with observations is at best highly 
speculative. 



4. Comparison with other work 

A detailed comparison between various evolutionary 
tracks available in the literature has already been done by 
Siess, Dufour & Forestini (2000). In the present section, 
we compare the evolutionary tracks most widely used by 
the community to describe the observational properties 
of objects in young clusters (Burrows et al. 1997, B97; 
D'Antona & Mazzitelli 1994, 1997, DM94, DM97; Palla & 
Stahler 1999, PS99). The comparison between these dif- 
ferent models in a theoretical HR diagram is illustrated in 
Figure |^ for a few low-mass star and BD masses. 

Let us first summarize the main differences between 
these models. One crucial difference is due to the outer 
boundary condition: DM94, DM97 and PS99 use approx- 
imate boundary conditions based on T(t) relationships 
assuming gray approximation and radiative equilibrium. 
Outer boundary conditions based on such approximations 
are wrong as soon as molecules form near the photosphere 
and convection reaches optically thin layers, i.e below 
T cff ~ 4000K (see Chabrier and Baraffe 1997, 2000; and 
references therein) They usually yield hotter effective tem- 
peratures for a given mass. Above ~ 4000^, the choice 
of the outer boundary condition is less consequential, as 
shown by the good agreement between the 0.8 Mq tracks 
of BCAH98 and of PS99 (see Figure |), all the other physi- 
cal inputs being similar for such mass. Note that, although 
the PS99 models start from a more realistic birthline, this 
does not affect the results after 1 Myr, as expected from 
the tests on initial conditions performed in §3.1.1. 

The B97 models for low mass stars and hot brown 
dwarfs (T G ff > 2000K) are based on gray atmosphere mod- 
els obtained by solving the radiative transfer equation, as 
described in Burrows et al. (1993). Such an approximation, 
although it represents an improvement over the previous 
T(t) relationships, still overestimates T e ff at a given mass 
compared to models based on full non-gray atmosphere 
models. This is certainly the reason why for VLMS on the 
Main Sequence (for m > O.O8M ), the B97 models are 
about 100-200K hotter than the BCAH98 ones. We note, 
however, that the very early evolution of the B97 mod- 
els proceeds at much hotter T c g than any other model. 
Although also true for masses below 0.2 Mq, this does 
not appear in Fig. || since the tracks start at 1 Myr and 
the 0.06 Mq of B97 displayed in Fig. || is indeed much 
hotter but for t < 1 Myr. The reasons for such differences 
along the Hayashi line is not clear. We only notice that 
similar shape of the Hayashi line (i.e. a strong decrease 
of T c g) is obtained if one assumes a fully adiabatic initial 
structure and adiabatic convection (i.e a m i x — > oo). Based 
on test cases, we find that such initial models are much 
hotter for a given L. If T Q g is high enough for H2 forma- 
tion to be negligible, such models start to evolve toward 
cooler T c ff. Once the fraction of H2 becomes significant 
and the adiabatic gradient small enough near the photo- 
sphere, evolution proceeds almost vertically in the HRD, 
as also expected from the analysis of Hayashi & Nakano 
(1963) (see §3.1.1). Such an evolution resembles the shape 
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of the B97 tracks. It may however be a pure coincidence, 
since B97 do not describe their initial conditions. 

Convection in all afore-mentioned models is treated 
within the framework of the mixing length theory, with 
1 < «mix < 2, except DM97, who use the Canuto- 
Mazzitelli formalism. As demonstrated in §3, the treat- 
ment of convection affects significantly evolutionary tracks 
and a variation of a m i x modifies the shape of the Hayashi 
lines (as long as super-adiabatic layers are present). Figure 
^ displays quite similar shapes for a given mass between 
BCAH98, PS99 and DM94, although not at the same T cS 
because of different outer boundary conditions. The DM97 
Hayashi lines behave differently with respect to models of 
other groups, as already noticed by Siess et al. (2000). 
This is certainly due to their different treatment of con- 
vection. Unexpectedly, the DM97 models predict a MS for 
VLMS close to the BCAH98 MS, which, as mentioned by 
DM97, is purely coincidental and stems from unexpected 
canceling effects between different treatments of convec- 
tion and outer boundary conditions. This is illustrated for 
the 0.2 Mq in Figure g. Since the grey-like outer bound- 
ary condition used in DM97 is expected to yield hotter ef- 
fective temperatures than the BCAH98 models, the close 
agreement on the MS suggests that the Canuto-Mazzitelli 
convection treatment yields larger super-adiabatic layers 
in the atmosphere. This is indeed required to decrease T e ff 
for a given mass. Such a behaviour, however, is not found 
by the recent simulations by Ludwig et al. (2001) (see 
§3.2). Indeed, the Canuto-Mazzitelli treatment of convec- 
tion yields results at odds with 3D hydro simulations for 
the outer thermal profile of the Sun (see e.g. Nordlund & 
Stein 1999), and does not provide an accurate treatment 
of convection in optically-thin media, at least for solar- 
type stars and low-mass stars ( see e.g. §2.1.3 of Chabrier 
& Baraffe, 2000 for a discussion of this topic). 

5. Observational tests 

In the previous sections, we emphasize the large uncer- 
tainty of evolutionary models at early ages and low grav- 
ities. Testing different initial conditions, we find out that 
after a few Myr these uncertainties become inconsequen- 
tial and all our models converge toward the same posi- 
tion in a HRD at a given age, for a given mass. Such a 
result does not necessarily imply that the models are reli- 
able at ages of a few Myr, since we have only tested sim- 
ple cases. More sophisticated initial conditions are beyond 
the scope of the present paper and were already explored 
by Hartman et al. (1997), illustrating the sensitivity of 
tracks and birthline positions to the (poorly known) de- 
tails of the protostellar accretion process (geometry, rate, 
temperature of added matter). 

A better knowledge of initial conditions may come 
from the determination of the minimum age below which 
present models start to depart significantly from obser- 
vations. Estimation of this age can constrain the charac- 
teristic time-scales and accretion rates of the protostel- 
lar collapse phase. Unfortunately, direct comparisons of 




Fig. 8. Comparison of evolutionary tracks from different 
authors, as indicated on the figure. All models start at t 
= 1 Myr, and evolve up to the MS for the 0.8 and 0.2 M Q 
stars. 

observations with models directly in colour - magnitude 
diagrams are extremely uncertain due to the large extinc- 
tion in star formation regions, which affects the observed 
energy distribution and thus the spectra and the colors. 
Only very few exceptions, such as a Orionis, exhibit low 
extinction. Recently, Bejar et al. (1999) and Zapatero et al. 
(1999, 2000) obtained optical and near-IR photometry for 
low mass objects in this cluster. In a (/ — J) vs M\ CMD, 
the data lie between the 1 and 10 Myr isochrones, respec- 
tively, for masses down to ~ 0.01Af Q , using the BCAH98 
and CBAH00 models (Zapatero et al. 2000; Bejar et al. 
2001). If statistics is improved and if the membership of 
the objects to the cluster is confirmed, such observations 
provide an unique opportunity to test directly the validity 
of young theoretical isochrones. They also offer the best 
chance to determine the Initial Mass Function (IMF) down 
to the substellar regime and the minimum mass formed by 
a collapse process (see Bejar et al. 2001). 

Young multiple systems provide also excellent tests for 
PMS models at young ages, because of the assumed co- 
evality of their different components. In addition, another 
strong constraint is supplied by the estimate of dynami- 
cal masses deduced either from binary systems (Covino et 
al. 2000; Steffen et al. 2001) or determined from the or- 
bital motion of circumstellar/circumbinary disks (Simon, 
Dutrey & Guilloteau 2000). An example is provided by the 
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quadruple system GG TAU (White ct al. 1999), with com- 
ponents covering the whole mass-range of VLMS and BDs 
from 1 Mq to ~ 0.02 Mq. Orbital velocity measurements 
of the circumbinary disk surrounding the two most mas- 
sive components imply a constraint on their combined stel- 
lar mass (Dutrey, Guilloteau & Simon 1994; Guilloteau, 
Dutrey & Simon 1999). This mass constraint and the hy- 
pothesis of coevality provides a stringent test for PMS 
models. The BCAH98 models are the most consistent with 
GG Tau (for details see White et al 1999; Luhman 1999) 
and provide the closest agreement with derived masses of 
other young systems (see Figj9j) . 

Most of the observed systems displayed in Figure [| arc 
better reproduced by tracks using a large value of a m i x 
(=1.9). However, for some systems, such as 1 (Covino et 
al. 2000), 2 (Steffen et al. 2001) and 4 (BP Tau from Simon 
et al. 2000), a better agreement is obtained with a m i X = 
1. Although a variation of a m i x with effective temperature 
and gravity is possible, as suggested by the simulations of 
Ludwig et al. (1999), none of these three systems occupies 
a peculiar position in (T c g, g) to suggest a different value 
of a m i x . This puzzle may reflect the uncertainties of PMS 
models based on arbitrary initial conditions. It may also 
be due to the large uncertainties of observationally-derived 
spectral type classifications, luminosity estimates and T e ff 
calibrations for such very young objects. Since they dis- 
play spectral features between that of giants and dwarfs 
(see Luhman 1999), a better representation of their spec- 
tral properties may require new indices more appropriate 
to these intermediate surface gravities, in the same vein as 
the pseudo-continuum ratios used by Martm et al. (1996) 
for Pleiades objects. The transformation of the inferred 
spectral type into T c g is even more difficult, because of 
the lack of reliable T c g - scales for such young T-Tauri like 
objects. Significant efforts were devoted within the past re- 
cent years to the elaboration of improved T c g - scales for 
M-dwarfs (Leggett et al. 1996) and M-giants (Perrin et 
al. 1998; van Belle et al. 1999). However, work remains to 
be done for T-Tauri like objects. Luhman (1999) defined 
a T e ff - scale intermediate between giants and dwarfs and 
based on the isochrone of BCAH98 which goes through 
the four components of GG Tau. Interestingly enough, 
applying this T e g - scale to young clusters such as IC348 
(Luhman 1999) and star forming regions like Chamaelcon 
I (Comeron, Neuhauser & Kaas 2000), the cluster mem- 
bers show a small scatter in age and no obvious correlation 
between age and mass. As mentioned by Comeron et al. 
(2000), this suggests in Chamaeleon I an almost coeval 
population which formed within less than 1 Myr. This 
short timescale supports the suggestion that star forma- 
tion is controlled primarily by large-scale turbulent flows 
rather than by magnetic processes such as ambipolar dif- 
fusion (Hartmann 2000). 

Given the assumption of coevality and the uncertain- 
ties of our initial conditions, the T e g - scale for young ob- 
jects suggested by Luhman (1999) needs to be confirmed 
in order to confirm these exciting results about star forma- 
tion process. Although still very preliminary, the compari- 



son of observed and synthetic spectra, as recently done by 
Lucas et al. (2001) for Orion objects, provides a promising 
way to define such a T e g - scale. 

6. Conclusion 

The very good agreement of models based on improved 
physics with observations for relatively old (t ;> 100 Myr) 
low-mass objects yields confidence in the underlying the- 
ory. Such evolutionary models can now be confronted to 
the complex realm of very young objects, providing impor- 
tant information on star formation processes and initial 
conditions for PMS models. Although based on extremely 
simple initial conditions (no accretion phase, no account 
of protostellar collapse phase and time scale, spherical 
symmetry), these models provide the most accurate com- 
parison with present observations of very young objects 
(dynamical masses, tests of coevality in multiple systems, 
CMDs). Given the combining effects of large observational 
and theoretical uncertainties at very young ages, however, 
one must remain cautious. It is probably too premature 
to conclude on the validity of the present models at early 
phases of evolution. We have examined in the present 
paper the uncertainties on evolutionary models of very- 
young low-mass objects arising from initial conditions, in 
particular the initial radius of the object, the efficiency of 
convection in the outermost layers and the initial abun- 
dance of deuterium. We have shown that at least the two 
first afore-mentioned uncertainties can affect drastically 
the fundamental properties, luminosity and effective tem- 
perature of objects younger than about 1 Myr. Therefore, 
any attempt to infer an age or a mass from observable 
quantities for these objects, in particular the initial mass 
function of very young clusters, must be considered with 
highly limited - if any - validity ! 

Realistic initial conditions can only be provided by 
multi-dimensional protostar collapse simulations, not by 
spherically-symmetric models for PMS initial conditions 
involving many free, ill- or unconstrained parameters. 
Because of numerical subtleties and complex physical pro- 
cesses (accretion fronts, turbulent time-dependent con- 
vection, hydrodynamical radiative transfer, magnetic field 
etc.), the construction of star formation models is a harsh 
task, which very likely will necessitate several years of ef- 
forts. Besides these theoretical difficulties, observations of 
very young objects can provide only limited guidance to 
such simulations, since most phases involved during the 
collapse are embedded in dusty cocoons. Only the final 
product can be observationally tested, when the proto- 
star becomes visible. This stage marks essentially the be- 
ginning of PMS evolutionary tracks. PMS tracks tested 
against observations thus provide a precious link to gather 
insight about star formation models from subsequent evo- 
lution. 

To progress in the field and in parallel with the devel- 
opment of star formation models, efforts can be directed 
toward: (1) a reduction of the main theoretical uncertain- 
ties affecting PMS models at low gravities, which involves 
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Fig. 9. Comparison of evolutionary tracks with observed PMS objects with derived masses. The BCAH98 tracks are 
displayed for 1.4, 1.2, 1, 0.8, 0.6 and 0.5 Mq (from left to right) for two values of a m i K . Observations are from Covino 
et al. (2000, filled triangles), Simon et al. (2000, crosses and open squares) and Steffen et al. (2001, open triangles). 
Open squares are single objects whereas all the other symbols indicate binaries. For each systems, masses are indicated 
(in Mq) on the right hand side of the figure. For GG tau A (9a+9b, crosses) the total mass is indicated. 




in particular a better determination of ^ m i x through multi- 
D hydrodynamical simulations and (2) the elaboration of 
a reliable T e ff-scale for young objects. 

Note: Tracks and isochrones for t > 1 Myr of the 
BCAH98 models (from 0.02 M to 1.4 Mq) and of the 
CBAH00 models including dusty atmospheres (from 0.001 
Mq to 0.1 Mq) are available by anonymous ftp: 

ftp ftp.ens-lyon.fr 

username: anonymous 

ftp > cd /pub/users/CRAL/ibaraffe 

ftp > get README 

ftp > get BCAH98_models.* 



ftp > get BCAH98_iso.* 
ftp > get DUSTY00_models 
ftp > quit 
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